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Turbulence-induced persistence in laser beam wandering 
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We have experimentally confirmed the presence of long-memory correlations in the wandering of a thin 
Gaussian laser beam over a screen after propagating through a turbulent medium. A laboratory-controlled 
experiment was conducted in which coordinate fluctuations of the laser beam were recorded at a sufficiently 
high sampling rate for a wide range of turbulent conditions. Horizontal and vertical displacements of the laser 
beam centroid were subsequently analyzed by implementing detrended fluctuation analysis. This is a very 
well-known and widely used methodology to unveil memory effects from time series. Results obtained from this 
experimental analysis allow us to confirm that both coordinates behave as highly persistent signals for strong 
turbulent intensities. This finding is relevant for a better comprehension and modeling of the turbulence effects 
in free-space optical communication systems and other applications related to propagation of optical signals in 
the atmosphere. @ 2015 Optical Society of America 
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The centroid of a laser beam experiences multiple deflec¬ 
tions when propagating through the Earth’s turbulent 
atmosphere due to stochastic refractive-index fluctua¬ 
tions along the optical path. As a consequence of phase 
changes, due to turbulent eddies with dimensions larger 
than the beam diameter, the laser suffers displacements 
perpendicular to the original direction of propagation. 
This phenomenon is known as beam wandering or spot 
dancing. In a first approach, the wandering of the laser 
beam in the turbulent atmosphere could be considered as 
fully random [1]. This means that the coordinate fluctu¬ 
ations of the laser centroid at a time are independent of 
those associated with previous instants. However, there 
exists some partial evidence that long-range correlations 
are present in the underlying temporal dynamics [2-8]. 
It is clear that a better understanding of the fluctuat¬ 
ing temporal behavior is essential for improving laser 
practical applications such as those related to tracking 
and communication purposes. Particularly, beam wan¬ 
der is considered the main cause of substantial signal 
losses at the receiver plane in free-space laser communi¬ 
cation systems, degrading its data transmission quality 
and reliability, and limiting its performance [9]. Indeed, 
beam-wander mitigating control systems have been re¬ 
cently proposed for overcoming this drawback [10,11]. 
Trying to shed some light on the laser beam wandering 
dynamics, in this Letter, we carefully analyze the tempo¬ 
ral correlations in the recorded position of a laser beam 
after it propagates through an indoor laboratory atmo¬ 


spheric chamber. We conjecture that this laboratory¬ 
generated turbulence is representative of fully developed 
atmospheric turbulence. Consequently, it is able to em¬ 
ulate the main properties of the turbulent flows that af¬ 
fect the laser beam of a communication link. In fact, 
similar devices have been used for testing atmospheric 
turbulence effects in several contexts [9,12-17]. Horizon¬ 
tal and vertical components of the centroid position of 
the laser spot are measured as function of time with a 
position sensitive detector located at the end of the prop¬ 
agation path. With the aim to statistically characterize 
the temporal correlations of these data streams, a de¬ 
trended fluctuation analysis (DFA) is implemented. DFA 
is a robust technique for detecting dependence among 
samples in noisy non-stationary time series. The results 
obtained from this fractal analysis confirm that the dis¬ 
placements of the centroid of the laser are consistent 
with long-memory correlated stochastic dynamics for the 
stronger turbulent conditions, i.e. when the turbulence 
effects can be resolved by the detector. 

A conceptually simple experiment was performed in 
controlled conditions in which a laser beam propagates 
through artificial turbulence—please see Fig. 1 for a 
schematic view of the optical setup. The wandering of 
the laser beam (10 mW HeNe Melles Griot Model 05- 
LHP-991) is detected by a position sensitive detector 
with an area of 1 cm^ (UDT SC-10 D). This detector 
measures the position of the centroid of the impinging 
laser beam with an accuracy of 2.5 /rm. Horizontal and 
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vertical coordinates were recorded at 500 Hz—we have 
confirmed similar findings for higher sampling rates. For 
the purpose of having fully developed inertial turbulence 
at stable and statistically repeatable conditions, we em¬ 
ploy a laboratory air turbulence generator, commonly 
called turbulator, similar to the one originally proposed 
by Fuchs et al. [18], and later enhanced by Keskin et 
al. [19]. To simulate the atmospheric turbulence, two air 
fluxes at different temperatures are forced to collide in 
the chamber producing an isotropic mix between hot and 
cold air. The hot source is an electric heater controlled 
by changing the current passing through it. The thin 
laser beam propagates across almost 0.35 m of turbu¬ 
lence in the mixing chamber. Air flow velocity is fixed 
because both fans operate at identical velocities so the 
turbulence characteristics are only due to the tempera¬ 
ture difference. By increasing the temperature of the hot 
source different turbulent intensities can be produced. 
The indoor chamber offers the advantages of full sys¬ 
tem characterization and repeatability in a single tur¬ 
bulent layer. The strength of the artificial turbulence, 
quantified through the structure constant C^, and the 
inner and outer scale were previously estimated following 
the procedure suggested by Masciadri and Vernin [12]. 
For such purpose, the variance of angle-of-arrival fluc¬ 
tuations of collimated laser beams as a function of the 
radius of different pupil masks passing through the tur¬ 
bulent layer is measured [12, Fig. 1]. By analyzing the 
averaging effect of the pupil sizes, and considering the 
effects of the inner {Iq) and outer (Lq) scales through 
the von Karman spectrum, it is feasible to estimate the 
different turbulence parameters (C^, Iq and Lq) by the 
fit of the theoretical model [12, Eq. (12)] to the empir¬ 
ical variances. Finally, is expressed as a function of 
the temperature difference between hot and cold sources 
(Ti and T 2 , respectively, in Fig. 1). Experiments with 
twelve temperature differences AT = Ti — T 2 ranging 
from 5 to 180 °C were carried out. Since turbulence is 
based on temperature gradient and air flow mixing, ref¬ 
erence measurements were taken in two different condi¬ 
tions: fans on and fans off, both with the heater discon¬ 
nected. In particular, measurements with the fans off can 
be considered as a background measurement that quan¬ 
tifies the electronic noise and room turbulence effects. It 
is worth noting here that the estimated structure con¬ 
stants are, at least, two or three orders of magnitude 
larger that those expected in outdoor experiments. Since 
the small turbulent path length in the mixing chamber, 
these higher turbulence strength are required to become 
detectable the laser centroid fluctuations. 

DEA was introduced more than twenty years ago to 
characterize the fractal dynamics of a system from which 
a time series has been measured [20]. Being the most 
popular approach to detect the presence of long-term 
memory in data [21], at present there are more than 
1800 articles published on DEA and its applications 
according to the information extracted from the Sco¬ 
pus bibliographic database (accessed in March, 2015). 



Eig. 1. Schematic diagram of the laboratory experimen¬ 
tal setup. 

However, specific applications in the optical field are 
scarce [8, 22-24]. Briefly explained, the DEA method 
consist of five steps [25]. Eirst, given a time series 
S = {xt^t = 1,...,A}, with N being the num¬ 
ber of equidistant observations, the cumulated data se¬ 
ries Y ii) = Y^\^i{xt — {x)), with i = 1,...,A and 

(^) = (Sfci /^■> is considered. In the second step, 
this profile is divided into [A^/sJ nonoverlapping win¬ 
dows of equal length s ([aj denotes the largest inte¬ 
ger less than or equal to a). A local polynomial fit 
Uv, m(i) of degree m is fitted to the profile for each 
window V = I,...,[A/sJ as the third step. The de¬ 
gree of the polynomial can be varied to eliminate con¬ 
stant (m = 0), linear (m = I), quadratic (m = 2) or 
higher order trends of the profile. Then, in the fourth 
step, the variance of the detrended time series is evalu¬ 
ated by averaging over all data points i in each segment 

= (I/s)ELi {F[{.v+ 

for V = !,...,[ A/sJ. In the last step, the 

DEA fluctuation function is obtained by averag¬ 
ing over all segments and taking the square root, 

Fm (s) = |(I/ LA^/sJ) K {v, s)] } ^ . This pro¬ 

cedure should be repeated for different values of the time 
scale s in order to unveil the s-dependence of Fm- If the 
time series has long-range power-law correlations, Fm (s) 
scales as 

Fm (s) (1) 

for a certain range of s. The Hurst exponent, i.e. the 
scaling exponent iJ, is estimated by the slope of the best 
linear regression in a double logarithmic plot. It quan- 
tihes the long-range correlations embedded in the time 
series: when H > 1/2, consecutive increments tend to 
have the same sign so that these processes are persistent. 
For iJ < 1/2, on the other hand, consecutive increments 
are more likely to have opposite signs, and it is said that 
the processes are anti-persistent. H = 1/2 is obtained 
for uncorrelated data [26, Chap. 9]. 

Twenty one independent realizations of 5000 coordi¬ 
nate points were obtained for each turbulent condition. 
Afterwards, the DFA analysis was performed for both 
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horizontal and vertical coordinates. Figure 2 shows the 
fluctuation functions obtained by implementing a DFA 
analysis with a detrending polynomial of second order 
(m = 2) for the different turbulent strengths. Represen¬ 
tative results for one particular realization for each tur¬ 
bulent condition are depicted, with increasing from 
bottom to top. A well-defined power-law behavior is con¬ 
cluded with higher slopes for the stronger turbulent in¬ 
tensities. Mean and standard deviation (over the twenty 
one independent realizations) of the Hurst exponents es¬ 
timated in the range s £ [30,1000] (vertical dashed lines 
in Fig. 2) as a function of for both coordinates are 
plotted in Fig. 3. In order to better interpret the re¬ 
sults obtained for the Hurst exponent. Fig. 4 shows the 
signal-to-noise ratio (SNR) for the different turbulent 
conditions. SNR was estimated as the variance of the 
signals at turbulent states relative to the value associ¬ 
ated with the background measurements. From Fig. 4 it 
is concluded that higher are needed to fully resolve 
the turbulence effects in the turbulator with the imple¬ 
mented detection system. 

For the reference measurements the detector is un¬ 
able to resolve the position differences, and a fully un¬ 
correlated electronic noise associated with the detector 
is measured. Consequently, the Hurst exponent is near 
0.5 as expected. As the turbulent strength increases, the 
SNR is larger than one (please see Fig. 4) and the de¬ 
tector begins to discern the turbulence influence. Simul¬ 
taneously, the Hurst exponent shows an increasing be¬ 
havior, for both horizontal and vertical coordinates, sat¬ 
urating at a value close to 5/6 for the higher values 
(please see Fig. 3). Results obtained allow us to confirm 
that turbulence introduces memory effects in time se¬ 
ries wandering because highly persistent dynamics are 
clearly concluded for the stronger turbulence intensities. 
Moreover, the similarity between horizontal and vertical 
estimated Hurst exponents confirms the isotropy of the 
turbulence within the laboratory chamber. It is worth 
mentioning here that H = 5/6 has been originally pro¬ 
posed theoretically for the turbulence-degraded wave- 
front phase within a Kolmogorov model [27], and very re¬ 
cently confirmed experimentally for the angle-of-arrival 
fluctuations of stellar wavefronts propagating through 
atmospheric turbulence [24]. 

Summarizing, a persistent stochastic fractal behav¬ 
ior is clearly concluded from the DFA analysis of laser 
beam wandering in laboratory-generated turbulence. Es¬ 
timated Hurst exponents, for both coordinates, converge 
to a value close to 5/6—the theoretical value associated 
with Kolmogorov turbulence—for the stronger turbulent 
conditions. Outdoor similar wandering measurements for 
horizontal paths in different moments of the day are be¬ 
ing planned for the near future in order to confirm the 
presence of memory effects in real atmospheric channels. 
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(UNLP), Argentina; Comision Nacional de Investigacion 




Fig. 2. Fluctuation functions ^ 2 ( 5 ) as a function of the 
scale s for the horizontal (upper plot) and vertical (lower 
plot) coordinate fluctuations for the different turbulent 
conditions. Results obtained for one particular realiza¬ 
tion of the twenty one recorded are depicted. A detrend¬ 
ing polynomial of order m = 2 and 96 different scales 
s £ [10,Af/4] equally spaced in the logarithmic scale 
were employed in the DFA implementation. increases 
from bottom to top. The slope of the best linear fit 
obtained for each one of these fluctuation functions is 
the Hurst exponent estimator. Vertical dashed lines in¬ 
dicate the range in which the linear fit for estimating 
the Hurst exponent is performed. Straight (gray dashed) 
lines with slopes 1/2 (bottom line) and 5/6 (top line) 
are also shown as references. The behavior observed is 
representative for the whole data set and similar results 
are obtained for other detrendings (m = 1 and m = 3). 
The significant similarity between horizontal and verti¬ 
cal fluctuation functions can be attributed to isotropy of 
the artificially generated turbulence. 

Cientifica y Tecnologica (CONICyT, Chile) (FONDE- 
CYT 1140917) and Pontificia Universidad Catolica de 
Valparaiso (PUCV, Chile) (123.731/2014). 
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